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0^ , Abstract. We study the ergodic properties of quantized ergodic 

maps of the torus. It is known that these satisfy quantum ergodic- 
ity: For almost all eigenstates, the expectation values of quantum 
observables converge to the classical phase-space average with re- 
spect to Liouville measure of the corresponding classical observ- 
able. 

The possible existence of any exceptional subsequences of eigen- 
states is an important issue, which until now was unresolved in any 
example. The absence of exceptional subsequences is referred to as 
quantum unique ergodicity (QUE). We present the first examples 
of maps which satisfy QUE: Irrational skew translations of the two- 
torus, the parabolic analogues of Arnold's cat maps. These maps 
are classically uniquely ergodic and not mixing. A crucial step is 
to find a quantization recipe which respects the quantum-classical 
Q\ , correspondence principle. 

*^ In addition to proving QUE for these maps, we also give re- 

sults on the rate of convergence to the phase-space average. We 
give upper bounds which we show are optimal. We construct spe- 
cial examples of these maps for which the rate of convergence is 
arbitrarily slow. 
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^ ■ 1. Introduction 

H : 

One of the central problems of "Quantum Chaos" is the question of 
the asymptotic behaviour of eigenmodes of classically chaotic systems 
in the semiclassical limit. In particular, one wants to find their limiting 
"mass distribution" in a suitable sense. 

Consider for instance the geodesic flow on a compact Riemannian 
manifold M (or rather, on its co-tangent bundle), whose quantum 
Hamiltonian is, in suitable units, represented by —A, the positive 
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Laplacian on M. Let ipj be a sequence of normalized eigenfunctions: 
Atpj + Xjipj — 0, J M \ipj\ 2 — 1. A suitable quantity for measuring the 
concentration properties of the eigenmodes ipj, in both the position and 
momentum representations, is the distribution on the unit co-tangent 
bundle S*M given by 

(1.1) /6CTM)h(0 P (/%^). 

Here Op(/) is a zero-order pseudo-differential operator with principal 
symbol / G C°°(S*M) and Op is some choice of quantization from sym- 
bols to pseudo-differential operators. The operator Op(/) is a quantiza- 
tion of the classical observable /, and {Op(f)ipj, ipj) are the expectation 
values for the operator in the state ipj. 

A celebrated result in this direction is "Schnirelman's theorem''^ 
which says that if the flow is ergodic then these expectations converge 
to the phase-space average of the classical observable /, for all but pos- 
sibly a zero-density subsequence of eigenfunctions. This phenomenon 
is commonly referred to as quantum ergodicity^. There are no examples 
where it is known if there are any exceptional subsequences. The case 
where there are none is referred to as quantum unique ergodicity (QUE) 
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In this paper, we consider a compact model of the above situation, 
where the dynamics, instead of taking place in the co-tangent bundle, 
occurs in a compact symplectic manifold, namely the 2-torus T 2 . The 
(classical) evolution is then given by iterating a symplectic map of the 
torus. 

In order to quantize such a map, one looks for a Hilbert space of state- 
vectors of the system, which are required to be periodic in both position 
and momentum representations. This constrains Planck's constant h 
to be an inverse integer: h = 1/N, and then the state space Hn is 
finite dimensional, of dimension precisely N. The semiclassical limit 
is now iV — > oo. By means of an analogue of Weyl quantization, one 
defines quantum observables Op N (f) corresponding to smooth classical 
observables / G C°°(T 2 ). 

Given a symplectic map A of T 2 , the quantum evolution is given by 
specifying a unitary operator Un on the state space Hn, which satisfies 
a version of the "correspondence principle" (Egorov's theorem): 

(1.2) \\UN 1 Op N (J)U N -Op N {foA)\\ N -=??0, V/gC°°(T 2 ), 



1 announced in |L6| with full proofs given by Zelditch JlTj for hyperbolic surfaces 
and Colin de Verdiere Q in general, see also ||. 

2 There are other notions of ergodicity in quantum mechanics, such as von Neu- 
mann's ifhil n2|, which are not related to the one used here. 



QUANTUM UNIQUE ERGODICITY FOR PARABOLIC MAPS 



3 



where foA{\) = f(A(q)), that is one requires that in the semiclassical 
limit, quantum evolution becomes classical evolution. The analogue of 
eigenmodes are then the eigenfunctions of the propagator Un- 

The main focus in the literature has so far been on hyperbolic trans- 
formations of the torus, the so-called cat maps |7], [TTJ, [|, to which 
the proof of Schnirelman's theorem [T6|, [T7|, [TJ can be adapted to prove 



quantum ergodicity, but not QUE [fj, [LB]. Assuming the Generalized 
Riemann Hypothesis, Degli Esposti, Graffi and Isola || found an ex- 
plicit infinite (though sparse) subsequence of values of N, for which 
they show that the expectation values for all eigenfunctions ip G Hn 
converge to the phase space average. 

Here, we will study a parabolic map of the torus (also called a skew 
translation), which is specified by choosing a real number a, and then 
defining^ 

a ° : (?) ~ G p + + 2 ;) - * 

When a = the motion is clearly integrable as p is a constant of 
the motion. For rational values of a, the map is "pseudo-integrable" 
in that the dynamics of the map on an orbit can be identified with an 
interval exchange transformation. For a irrational, the map is ergodic 
and in fact it was found by Furstenberg |§ to be uniquely ergodic. 
These maps possess no further "chaotic" properties; for instance they 
are not mixing. 

We propose a quantization procedure that at each value of N, re- 
places a by a rational approximant a/N, and then construct a unitary 
propagator lA a ,N on Ti^ which satisfies an exact version of Egorov's 
Theorem: 

U a^N Op7v(/)W a ,iV = Opjv(/ o A a/N ). 
Then taking any sequence a/N a we show that ( |1.2|) holds. This 
gives us a quantization of the map A. There are other recipes in the 
literature |3|, |2|]; however, they do not satisfy (|1.2|). 

Once we have the analogue of Egorov's theorem ( |1.2| ) and have set up 
the necessary tools from pseudo-differential calculus on T 2 , we show: 

Theorem 1.1 (QUE for parabolic maps). Suppose a is irrational, f G 
C°°(T 2 ) a smooth observable, and a/N — > a a sequence of rational ap- 
proximants. Then for any normalized eigenfunctions ip G TLn of the 
propagator Ua^N, we have 

(Op N (f)^p,ip) -> / f{p,q)dpdq, N -> oo. 



3 For a technical reason we shift q by 2p rather than p. 
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That is the parabolic map A a satisfies quantum unique ergodicity. 
This is the first known example of QUE. 

The remainder of our paper concerns the rate of convergence. We 
take approximants so that \a — a/N\ < 1/N. Suppose that a is badly 
approximable (in the sense that \a —p/q\ ^> e q~ 2 ~ e for all e > 0). We 
then show (Corollary |4.3| ) that for any normalized eigenfunction ip of 
the propagator U a ^ we have 

KO Pjv (/)^>- / f(p,q)dpdq\<^ u N~ 1 ^, Ve>0. 
Jt 2 

The reason why the rate is _/V _1 / 4+e and not, as one might have guessed 
are degeneracies j n the spectrum, which occur whenever a and 
N are not co-prime. We can, however, always construct an explicit 
basis of eigenf unctions ipj (j = 1, . . . , N), for which 

|(0 PiV (/)^,Vi}- / f(p,q)dpdq\<^ f N- 1 / 2 , 
Jt 2 

see Section |5|. In the absence of degeneracies (a and N co-prime), we 
thus indeed obtain a rate of N~ l l 2 (cf. also Theorem 4.1). 



As for lower bounds on the rate, we show (Theorem |6.1| ) that for 
the observable f(p, q) = e 2nv2p , for all irrationals there is a sequence of 
values of iV and normalized eigenfunctions ip G 7Yat for which 

1 



\(Op N (f)M)- f(p,q)dpdq\^> 



T 2 



ATI/4- 

Thus for badly approximable a, Corollary [O] is sharp. Moreover, 
unlike the situation with badly approximable a, we can construct irra- 
tionals for which the rate of convergence in Theorem is arbitrarily 
slow, e.g. slower then 1/ log log log iV (Theorem 
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2. Quantum mechanics on T 2 

2.1. Notation. We write e(x) = e 2mx and e^aO = e^ x . Z^v denotes 
the residue class ring Z/iVZ. A <C e B and A = O e (B) both mean that 
there is a positive constant c depending only on e, such that \A\ < c|S|. 

2.2. The Hilbert space of states. To recall the basics of quantum 
mechanics on the compact phase-space T 2 [7|, f|, we begin by de- 
scribing the Hilbert space of states of such a system. The guiding rule 
is Heisenberg's uncertainty principle, which asserts that simultaneous 
measurements of momentum p and position q of a quantum particle are 
ambiguous within Planck cells of volume h (Planck's constant). Hence 
if the phase space volume V is finite, the dimension N of the Hilbert 
space TCn describing the state of the quantum particle has to be finite 
as well, and is precisely given by N = V/ h. 

In the case of the torus T 2 , we take state vectors to be distribu- 
tions on the line which are periodic in both momentum and position 
representations: ij){q + 1) = i/j(q), [J'hi^ip + 1) — [^ 7 hi ) ] (p) > where 
[T h ip]{p) = h~ 1 ' 2 J i)(q)e(-pq/h)dq. The space of such distributions 
is finite dimensional, of dimension precisely iV = 1/h, and consists of 
periodic point-masses at the coordinates q = Q/N, Q G Z. 

We may then identify Ti^ with the iV-dimensional vector space 
L 2 (Zjv), with the inner product ( • , • ) defined by 

(2-1) M) = ^ £ <KQMQ), 

QmodiV 

This inner product induces a norm ||-|| on the space of operators on 
T~Cn, that is on the space of N x N matrices. 

The Fourier transform jFjy may now be defined as the unitary map 

(2.2) $(P) = [F N i(,](P) = N-i HQ)e N (-QP), 

QmodN 

its inverse F^ 1 is then clearly given by 

(2.3) 1>{Q) = tt x $\{Q) = N-* $(P)tN(PQ). 

PmodN 

2.3. Translation operators. A central role will be played by the 
translation operators 

M](Q)=V(Q + 1) 

and 

M](Q)=e N (QMQ), 



6 JENS MARKLOF AND ZEEV RUDNICK 

which may be viewed as the analogues of differentiation and multiplica- 
tion (respectively) operators in usual Fourier analysis on M. n . In fact in 
terms of the usual translation operators on the line qip(q) = qip{q) and 
pip(q) = -^[^^{q)i they are given by t\ = e(p), t 2 = e(g). Heisenberg's 
commutation relations read in this context 

(2.4) t a 1 t\ = t b 2 t a 1 e N (ab) \/a,beZ. 

The Fourier conjugates of t\ and t 2 are 

and 

T^t 2 T ^y 1 —t 1 1 . 
2.4. Observables. For n = (ni,n 2 ) E Z 2 put 

T N (n) = e N (^)t n 2 X\ 

Then 

T N {m) T N {n) = e N { y ^ ! ) T N {m + n) 

with the symplectic form 

uj{m, n) = m\n 2 — m 2 n\. 

For any smooth function / E C°°(T 2 ) on our phase space T 2 , define 
a quantum observable 

Op N (f) = Yl Kn)T N {n) 

where f(n) are the Fourier coefficients of /. The observable Op N (f) is 
also called the Weyl quantization of f. 

We have Op N (f)* = Op N (f) and hence Op N (f) is self-adjoint for 
real- valued /. 

The connection of these quantum observables with the "classical" 
translations of the torus 

is explained in the following lemma. 

Lemma 2.1. For every f E C°°(T 2 ) ; a E 7h, we have (i) 

tiO Pw (/)tr a = o Pjv (/o5 2 a/Ar ), 
t a 2 op N (f)t 2 a = op N (f°sr /N ), 
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and (ii) for all a G K, 

||t? Op N (fK a ~ Op N (f o S%)\\ « 7 |a - £|, 

¥2 Op N (f)t 2 a ~ Op N (f o 5r Q )|| « f | a - 



Proof. With the commutation relations (|2.4J) we find 
ttOpM)^ = Y,K n yN{an 2 )T N {n). 

n 

On the other hand, 

Op N (f ° S2) = E 7He(tm2)T JV (n), 

n 

and the bound 

|e(an 2 ) - e (^™2)| < |27rn 2 ||a - -j^| 

concludes the proof of the statements concerning £1. The results for £2 
follow accordingly. □ 

2.5. Priedrichs symmetrization. Let h G iS(1R 2 ) be an even, real- 
valued Schwartz function normalized such that 

h(x) 2 dx = 1. 

JR 2 

The kernel 

Kat(x, a/) = ^2 h ( N ^ (x-x' + m)) 

is now used to define an alternative quantization (a variant of the "anti- 
Wick quantization") 

Op N (f) ■ L 2 (M -> L 2 (Ziv) 
of the observable / G C°°(T 2 ), by setting 

Op N (f) = ^[ [Op N (K N (-,x))] 2 f(x)dx. 

Jf 2 

The normalization constant 



c N = ^2 / h ( x + rfx 

„^-77? JR 2 



n£Z 2 

is chosen such that 

(2.5) Qpjv(l) = idjv. 

Asymptotically, 

C N = 1 + R (N- R ) } any/?. 
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The main feature of this quantization is positivity: If / > then 

(2.6) {6p N (f)M)>0, 
since 

(0 VN (f)^^) = -L [ \\Op N (K N (.,x)W 2 f(x)dx, 
which is clearly non-negative. Hence 

defines a measure on C°°(T 2 ), with total mass 



This "positive" quantization differs from the Weyl quantization at 
most by terms of order 1/N, as stated in the following proposition. 

Proposition 2.2. For every f G C°°(T 2 ) we have 

\\0p N (f)-0p N (f)\\^f^- 

Proof. By the Poisson summation formula, our kernel can be re- 
expressed in the form 

K N (x, x 1 ) = — - V h{— -) e(m(x - x')), 

N2 ^ 2 iV2 

where h is the Fourier transform of h. Then, by definition, 
Op N (K N ( ■ , x)) = V h(-^r) e(-mx) T N {m) 

and 

[Op N (K N (.,x))} 2 = ±- M^)^) x 

iv ^— ^ TV 2 iv 2 



x e(— (m + n)x) e.w( ) T N (m + n) 



With this, we find 



AT ' V «2 \" ^ 



= £ /(A;)G w (- r )T J v(A;) 



with 

, s 1 \ -v m m , ,1 , m ss 

= -— J] ft*- e w t . 

TVCat ^— ' TV 2 TV 2 2 TV 2 
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Therefore 

Op N (f) ~ Op N (f) = E ( Gn ^ ~ Tn ^- 

We have Gat(O) = 1 by Poisson summation and the definition of Cat. 
It is easy to see that Gn and its derivatives are bounded uniformly in 
N by rapidly decreasing functions of t. Moreover, Gjv(— t) = Gjv(t) is 
even as is easy to see using h is even and the bilinearity of uj. Thus 
expanding Gjy(t) in a Taylor series at t — and noting that since Gn 
is even, the first order terms are missing, we find that for \t\ <C 1, 

G N (t) = G N (0) + O(\t\ 2 ) = l + O(\t\ 2 ). 
Therefore since the Fourier coefficients f(k) are rapidly decreasing, 

l|OPiv(/) - op N (f)\\ « E |7(fc)|| Gjv ( * ) - 1| 

« E E 

|fc|<JVVio ' |fe|>ArVio 



1 1 1 

and the Proposition follows. □ 
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3. Quantizing skew translations 

In this section we define a quantization U a ^ : L 2 (Zjv) — > L 2 (Z N ) for 
the skew translation of the torus 



A a : 



qj \q + 2pj' 



We define the quantization in the momentum representation, that is 
Ua,N = ^ r jv 1 V a ,AT^Ar, by choosing an approximation a/N to a, with 

a . 1 

(3-D 1-jfKjf 
and then setting 

[V a ,Nip](P) := e N (-(P - a) 2 ) ^{P - a). 

The relation between the quantized map and the classical map A a 
is given by 

Theorem 3.1 (Egorov's Theorem for A a ). If \ot — a/N\ < 1/N then 
for every f e C°°(T 2 ) we have 

\\K,n Op N (f)U a>N - Oj> N (f ° A a )\\ </-. 

This is an immediate conclusion of the following Proposition together 
with the choice (13. II): 



Proposition 3.2. For every f £ C°°(T 2 ) we have (i) 
(3-2) U-jf VN {f)U a , N = Op N (f o A a/N ), 

and (ii) for all a£M, 

\\K,n Op N (f)Ua,N-Op N (foA a )\\ |«-^| 

Proof. We have to show that 

Va,N OPiv(/) Va,N = 6p N (f ° A a/N ), 

and that for real a, 

HV-^Op^Va^-Op^/o^H «/ \a-jj\ 

where 

6p N (f)=^NOp N (f)^N- 

Note that we can write 

Va,N=ti a V0,N, A a = S?oA . 
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Since 

Kir Op N (f) V a>N = V^ N t a x F N O VN {f)^\ a V 0>N 

= Kn^ n ^ a QPNif^F^Vo^N, 

we find, by virtue of Lemma |2.1| , 

K,N OPat(/) V a , N = V^ N PN (f o S? /N ) V ,N, 

and 

\\KN^NWKN-Vo, l N 6v N (f°s?)VoM\ «/ 

respectively. It thus remains to be checked that 

v ^Op 7V (/)v ^ = Op JV (/oA ). 

To this end, note first that 

n 

Second, let us show that 

(3.3) trX 1 V , N = e N {nl) t~ n H n 2 ^ 

holds: 

= e N (P 2 )[t^V , N MP-^) 

= e N (P 2 + ni (P - n 2 ) - (P - n 2 ) 2 )^{P - n 2 ) 

= e N {n 2 2 )e N ({ni + 2n 2 )(P - n 2 ))ip(P - n 2 ) 

= e N (n 2 2 )[t? +2n ^)(P - n 2 ) 

= e N {nl)[r l n Hf +2n ^]{P). 

The commutation relations (|3.3|) now lead to 

Kn°pMKn = Y.f^ -2n 2 ,n 2 )e N (^)t; n H?. 

n 

The Fourier coefficients of / o A are, however, exactly /(ni — 2n 2 , n 2 ), 
and our proof is complete. □ 



Remark 3.1. There are other quantization schemes of skew translations 
in the literature Q . However, they do not satisfy Theorem |0| and 
so their relevance to the classical dynamics is unclear. 
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3.1. Proof of Theorem |1 . 1|. For each N, choose an approximant 



a/N with \a — a/N\ —> 0, and a normalized eigenfunction ip G 7~In of 
U a ,N- Using the results of Section [27|, we get a sequence of probability 
measures (J>N,ip{f) = (Op N (f)ip, ip). Since they differ from the distribu- 
tions / i— > (Op J y(/)'0, '0) by terms which vanish as — * oo, it suffices 
to show that converge to Lebesgue measure A : / i— > J T2 /. 

To see this, note that the space of probability measures on T 2 is com- 
pact, and hence any sequence of probability measures has a convergent 
subsequence. Thus the sequence /ijv,i/> has a limit point, which is a prob- 
ability measure. Any such limit point v is then invariant under the map 
A a by Egorov's theorem (Theorem |3.1| ). For irrational a, the map A a 
is uniquely ergodic which forces v — A. Thus Lebesgue measure A is 
the unique accumulation point of our sequence. This forces fijf,ip A, 
otherwise there would be a neighborhood of A which excludes infinitely 
many fj,^. But then these latter would have to contain a convergent 
subsequence whose limit would not be A — a contradiction. 
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4. Upper bounds for the rate of quantum unique 

ergodicity 



Besides the convergence result of Theorem O , we can also give a 
bound for the rate of convergence. To do this, we will always assume 
that we pick approximants such that \a — a/N\ < 1/N. Our first result 
is 

Theorem 4.1. If a is irrational then for all f G C°°(T 2 ) ; and any 
normalized eigenfunction ip G L 2 (Z N ) of the propagator lA a ^ , 

(Op N (m 1>) = f f{P, q) dpdq + O f (j^), iV ^ oo, 

where M = N/ gcd(a, N) . 

To see that this has content, we note 

Lemma 4.2. Suppose we take approximations a/N to a with — > a 
as N —>■ oo. If a is irrational then N/ gcd(a, N) —>■ oo as N — ► oo. 

Indeed, write D = gcd(a, N), M = N/D, b = a/D. If a is irrational 
then c M '■= min{|a — ^| : m < M} > 0, and so \a — fj \ = \a — jj\ > 
cm > 0. Thus if M is bounded then cm is bounded away from zero, 
contradicting £ — > a. 

We say that an irrational a is badly approximable if 

a , 1 
a --» e -5Z7> Ve>0. 
n n /+e 

In that case, we can say something stronger then just that M — > oo as 
N — > oo. In fact we have M 3> e jV 1//2_e since 

1 a . . b . 1 

— > lot — — \ — lot — — 



N N M M 2 + e ' 

We thus find: 

Corollary 4.3. // a is badly approximable then for any normalized 
eigenfunction ip of the propagator U a ,N we have 

(Op N (f)i;,i>) - [ f\p 1 q)dpdq<^ u N- 1 /^ foralle>0. 

4.1. Proof of Theorem ¥7T. The idea of the proof of Theorem |4.1| is 
to use the fact that for an eigenfunction if) ofU a ^, we have (Op N (f)ip, if)) - 
{Op N (f T )^^), where f T := ^Y^=xf ° K/n is the ergodic average 
of /. Taking T = N we show directly that for any ip G 7Yat, we have 

KOp^n^^-^/KM- 1 ^. 

We start the argument by taking for / the basic exponential e mtU (p, q) 
e(mp + nq). 
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Lemma 4.4. For any normalized ip 6 L 2 (Z N ), 
(4.1) 

i/r> / t \ / ,\i2^ 1 \7( \\i\Q( an a(m - n) + n 2 + 2ny 

\(Op N {e m ^,^)\ }^ my)\ \ s (-tf> jj ; T )I 

ymodN 



2 



w/iere 

T 

(4.2) S(a,b;T) :=^e{at 2 + bi). 



t=i 



Proof. Iterating A a gives 

(«) <f!) = f. i ' + S,J-dl. 



From (14.31) we find that 



(4.4) e m , n o A l a/N = e N (ant 2 + a(m - n)t))e m+2 tn,; 



%/N 

Therefore 

T „ T 



e m,n : = f Yl e ' m > n ° = f ^{ant 2 + a(m - n)t))e m+2t n,n, 



„T _ 
,n ' ^-i 

t=l " /-I 

and quantizing we get 

T 



(4.5) Pjv (e^J = i ^ e N {ant 2 + (m - n)t))T JV (m + 2tn, n) 



T 

t=i 



In particular, ||Op 7V (e^ „) || < 1. 
From ( [4.5|) we get 



— e^v f cmt 2 + a(m — n)t H — (m + 2tn)n + nxj ip(x + m + 2tn), 
t=i \ " ' 



>N\ cT 

T 

T 

and so 



_ / ran ~> 



xmodN 



c 7\r (ant 2 + (a(m — n) + n 2 )t) ip{x + m + 2tn). 
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On applying Cauchy-Schwarz we find 
\(Op N (e T m , n )^)\ 2 < ^ £ |V>(*)| 2 

xmodN 

1 T 

•— I'^^e^ (ant 2 + (a(m — n) + n 2 )t) tp(x + m + 2tn)\ 2 . 

xmodN t=l 

Now ±£ 

xmodN I ^0*0 I — IIV'II ~~ 1 an d using the Fourier expansion 
ip(x) = N~~ l l 2 ^2 v 4>(y)e(yx) in the second x-sum gives 



|(0 Pjv (e^ n )V^)i 2 
1 1 



x,y,y'modN 

r 



e N (an(t 2 - t' 2 ) + (a(m - n) + n 2 )(t - t') + 2nyt - 2ny't') 



t,t'=\ 

r 



= T^jV ^ \^{y)\ 2 ^^e N {ant 2 + (a(m~n)+n 2 + 2ny)t)\ 2 , 

ymodN t=l 

by Parseval's identity. □ 

We now take T = N and then get a Gauss sum for S(a, b; T) in 
Lemma O : Define the complete Gauss sum 

G(c,d;N) := £ e7V (ct 2 + dt). 

tmodTV 

We will need a very classical estimate of its absolute value, which we 
recall: 

Lemma 4.5. If gcd(2c, N) = 1 then 

\G(c,d;N)\ =N 1 ' 2 . 

If c 7^ mod iV then 

\G(c,d;N)\ < iV 1/2 gcd(2c,iV) 1/2 , c^O mod N. 
If c = mod iV i/ien 



G(0, d; N) 



N, d = mod N 
0, d ^ mod N. 
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Proof. Since the case c = mod N is obvious, we assume c 7^ mod N. 
By multiplying out \G(c, d; N)\ 2 , changing variables and switching the 
order of summation we find 

\G(c, d; N)\ 2 = ~ *2) + d(h - h)) 

tifymodN 

= e N(cy 2 + dy) ^2 e N(?cyt). 

ymodN tmodN 

The inner sum is either iV or 0, depending if 2cy = mod N or not. 
This gives 

\G(c,d;N)\ 2 = N J2 e N (cy 2 + dy). 

y:2cy=0modN 

If gcd(2c, N) = 1 then the only solution of 2cy = mod N is y = 
mod N so we get equality \G(c,d; N)\ 2 = N, while in general the 
number of solutions is gcd(2c, N) which gives the bound \G(c, d; N)\ 2 < 
Ngcd(2c,N). □ 

Lemma 4.6. For any normalized if) G L 2 (Zjy), and \m\, |n] < M, 
(m,n) 7^ (0,0), we have, i/n^O, 

(4-6) \{0 VN {e N m ^)\<\2n\^ 2 M-^ 2 , 

while if n = but m/fl mod M then 

(4-7) (O PiV « )^> = 0. 

Proof. From Lemma |4.4j we have 

l(Op^(e^,^)| 2 <^ E my)\ 2 \G(an,a(m-n)+n 2 + 2ny;N)\ 2 . 

ymodN 

Recall that a/N = b/M with b, M co-prime, D = gcd(a, N). We have 
an ^ mod N if and only if n 7^ mod M. Thus if n 7^ mod M then 
by Lemma (L5] 

|(0 P7V (e^,^)| 2 <iV- 2 £ |^)| 2 gcd(2an,iV) 

ymodN 

= M~ 1 gcd(2n, M)||^|| 2 
< \2n\M-\ 



since II^H = HVII = 1. 
If n = then by ([L4 



- - ^ eTv(amt) 



e m,0 ~~ e m,0 

tmodTV 
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which vanishes if am ^ mod N, equivalently if m / mod M. 
Thus Op JV (e^ )0 ) = if m ^ mod M. □ 

4.2. Conclusion of the proof. If ip is an eigenfunction o{lA a ^, then 

r 



t=\ 
1 T 



t=l 

By Egorov (p]), 

Z<JOp w (/)^ = Op JV (/o4 /jy ) 

and so 

(op N (m^) = (op N (f T )^ip). 

Then expanding / in a Fourier series / = J2( m n) /( m ; n ) e m,n and 
applying the ergodic average operator with T = iV we get 

f N = Yl f^Kn- 

(m,n)^(0,0) 

Therefore 

(Op N (f N )M)- f f(p,q)dpdq= V /(m,n)(0 Pjv (ei n )^^)- 

/T2 ^ 

^ 1 (m,n)^(0,0) 

Now we have ||0 P j V (e^ n )|| < 1 and so we truncate the sum above to 
frequencies \m\, \n\ < M with error at most 

]T |/(m,n)|«M- x 

\m\>M or \n\>M 
(m,n)^(0,0) 

since f(m } n) is rapidly decreasing. It is important to note that since 
a is irrational, we have M — > oo as iV — > oo which we assume. 



For the small frequencies we use Lemma 4.6 to find 



?(m,n)(0p w (O,Vi>< £ |/(m,n)||n|^ M -^ 



|m|,|n|<M |m|,|n|<M 
(m,n)^(0,0) n^O 



Thus we find that for normalized eigenfunctions ip we have 

((W/)V,V>- / f(P,l)dpdq < f M" 1 / 2 + M" x « ikT 1 / 2 . 



18 
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This concludes the proof of Theorem 4.1 
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5. Explicit eigenfunctions 

We begin by calculating the eigenvalues and a basis of eigenfunctions 
of the quantum map V a ,N, defined by 

[V a ,Nip](P) = e N (-(P - a) 2 ) V(P - a). 
The eigenvalue equation 

(5.1) V a , N ^ = e N (<P)^P 

yields the following simple recursion relation for the eigenfunction ip, 

(5.2) ^(p + a ) = e7V (_0_p2 ) ^ (p)> 

We can now construct N linearly independent solutions ipj of ( |5.1|) , 
j = 1, . . . , N, as follows. Let D = gcd(a, N) be the greatest common 
divisor of a and N. Put 

a N 

6 = —, M=— , 
D D 

and write furthermore 

j = v + Di, v^ihD], ie[Q,M-i]. 

For a given j G [1, A/], the pair (rj, I) is uniquely determined. 
Proposition 5.1. The functions 

[ s/De N (- riau* - ulD + a 2 y (M-i)(2 A f-i)-(,-i)(2,-i) ^ 

if P = T) + ua mod iV 



if P ^ rj mod .D 
sofee Equation ( 5. 1\ ) with eigenphases 

, . , „(M-1)(2M-1) 

(5.3) = /p_^ + r?a _ a 2l ^ ; 

o 

and form an orthonormal basis o/L 2 (Zjy). 

Proof. Put ip v ,i(ij) = VD- The recursion relation (|5.2[) then implies 
that at points of the form P = rj + ua (u = 1, 2, . . . ) reads 

4>n,i(P) = VDe N ( - u(j) v> i - ^(77 + ma) 2 ). 

m=0 

Since A 7 " divides Ma, and thus ip v ,i(v + Ma) = ip v> i(rj), the eigenphases 
are determined by 

JW-l 

e N (M(j) v j+ J2(v + ma) 2 ) = 1, 

m=0 



MP) = { 
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leaving an ambiguity mod -p, which permits to put 

j M-1 



M 

m=0 



.2 ,x ,(M-1)(2M-1) 



= lD-r] 2 - r/a(M - I) - a 1 



6 

Since ?7aM = mod iV, we drop this term. A straightforward ma- 
nipulation leads to the expression for tp^i as given in the proposition. 
Orthonormality follows from 



and 



(*l>r,,U , *l>rf,v) = -jy VV( P )VV,*'( P ) 
P mod AT 
' M-1 

_£ ejv (-„(J-J')I>) if 17 = 1/ 
if T] 7^ 77' 



5?2> ("<'-'>) = o otherwise . 

u=0 k. 



□ 



Corollary 5.2. The multiplicity m(<f>) of an eigenphase <fi is bounded 
by m((f>) <C D2t(D) <C e D^ +e , any e > 0, where t(D) is the number of 
divisors of D. 

Proof. For a given 0, we would like to count the number of solutions of 

(5.4) W-^ + n a-a^ M - 1 ^ M - 1 K^ AN . 

6 

This implies that — rf = <p mod D. In order to count the number 
of solutions of the latter equation, define 



1 if x = mod D 
if x ^ mod D. 



S D (x) 

Then 

p=i 

Since 



pmodD 
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we find that 

pmodD ^ P=l ' 

The sum in brackets is a classical Gauss sum, whose absolute value is 
bounded by gcd(L>, 2p) (L emma 14.51), and thus 

< D3r(D). 

For fixed 77, Equation (|5.4j ) determines / uniquely mod M. □ 
Let us put 

m 2 = —nib, mi = (n 2 — 2pn — p(p — l)a)b — (/ — f). 



Lemma 5.3. We have 

n — n' = ni mod £) 



E^Lp 1 e M (m 2 ^ 2 +mit/) 



(J, otherwise 
and in particular for n\ = mod M 

Jo i/m^OmodM 
II mi = mod M . 

Proof. We have 

(Ttv^Vv^W') = ^ E e ^(~^ + n 2^)Vv( p + n i)VV,r( p ) 

PmodJV 

1 v > ,nin 2 , , w 1 / 1 \~r / 1 ^ 

= — 2^ ejv( — — + n 2(r/ + z/a)) ^(77 + 1/0 + n i) ^P v ',Av + va), 

which is non-zero only if there is a p such that 

(5.5) ni = pa + n — rf mod iV, 

hence in particular n\ + n' — 77 = mod £>. Taking absolute values and 
using the explicit expressions for the eigenfunctions we obtain 

(5.6) |(Tjv(n)Vv>Vy,*')l 

M-l r f+jtt-1 



1 

M 



^ejv n 2 va-((j) rj! i-4> rj > j if)u- ^2 {n+maf +^2{n' +ma) 2 



u=0 L m=0 m=0 

By virtue of Relation (|5.5|) we have 

^— 1 J/+/.1— 1 

/JC 7 / + ma) 2 = (77 + ma — ti^) 2 mod N. 

m=0 m=fi 
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This formula allows us to simplify ( |5.6| ) to 

(5.7) \{T N (n)ip^ h i) v >^)\ 



1 

M 



M-l 
u=0 



n^va — (4> v j — 4> v ' ,v ) v — niav(v — 1) — n x (2/ia+ 2t] — rii ) v 

M-l 



1 

M 



u=0 



with 



m 2 = —nib, 77i\ = (712 — 2fi7] — fi(fi — l)a)b — (I — I'). 
In the case t\\ = mod M we have m 2 = mod M and thus 



|(lV(n)Vv,Vy,J')l = 



M-l 



2j e MK 



if mi ^ mod M 

1 if rrii = mod M. 



□ 



In the sequel a G Z will be chosen such that 

a 1 

~ iv' < iV 



holds. 



Proposition 5.4. Lei / G C°°(T 2 ), and assume a is diophantine. 
Then 



(0 PjV (/)Vv,Vv>= / fdx + O f , a (N-*). 
If f is a polynomial, then the above relation holds for all irrational a. 
Proof. Without loss of generality we assume J j2 f dx = 0. We have 

(Ov N (f)ip v ,h Vv) = ^f( n ){ T N{n)'4> v ,h^r 1 ,i)i 

where f(n) are the (rapidly decreasing) Fourier coefficients of /. 

Following Lemma |5.3| , we distinguish two cases. 

Case A: ni ^ mod M. Then m 2 = — abii = —bni ^ mod M and 
by Lemma [4.5| 

M-l 

2j e A / (m 2 ^ 2 + miz/) 
from which we obtain 



< v /Mgcd(M,2|n 1 |) < fliji, 
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Case B. n\ = mod M. Here 

M-l 



\{T N {n)ij) v>h il}^i)\ 



M 



^2 e M{n 2 bv) 



u=0 



if n 2 ^ mod M 

1 if n 2 = mod M. 



In summary, 

(5.8) |(0 PiV (/)^,^}|< v5Jir*£)|n 1 |*|/(n)| + ^ |/( 



n 



ni^O 
Dim 



£>|ni 
M\ni,ri2 



the first sum corresponds to Case A, and the second sum to Case B. 
Since / is smooth we can bound the first sum by 

V2M~* |ni|'|/(n)| </, fl M-^D~ R = N~^D~ R+ ^ 

m 
D 

for any R. Hence 



D\m 



V2M~? KI*I/MI </ iV ~^- 



A' 



ni^O 
D|ni 

The second sum in (|5.8| ) is bounded by 

£ |/(n)|< />fl AT 

M|ni,ri2 

In particular the sum is empty for N large enough, if / is a polynomial, 
because for a irrational, M grows with N (see Lemma |4.2| ). Thus we 
get the second part of the proposition. □ 
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6. Lower bounds 

We begin with a result which implies that our bound on the rate of 
convergence (Corollary |4.3| ) for badly approximable a is the optimal 
one: 

Theorem 6.1. For any irrational a, there are arbitrarily large N, ap- 
proximants \oc — a/N\ < 1/N and eigenf unctions if) ofl4 a> N so that 

1(^(2,0)^)1 = ^^. 

Unlike badly approximable a's, where we have an upper bound on 
the rate of convergence of (Corollary pj), we can construct 

irrationals for which the rate of convergence is arbitrarily slow, e.g. 
slower then 1 / log log log N: 

Theorem 6.2. Let g(x) be an increasing positive function. Then there 
is an irrational a such that there are arbitrarily large values of N for 
which there are normalized eigenf unctions if) ofU a ,N satisfying 

|<7V(2,o)v,V}| 

6.1. Constructing special eigenfunctions. In order to prove The- 
orems and we first use the results of Section |5| to construct 
special eigenfunctions if) for which the upper bound of Theorem |4.1| is 
optimal: 

Proposition 6.3. If D = gcd(a,iV) > 2 and M = N/gcd{a,N) is 
odd then there are normalized eigenfunctions if> so that 

(6-1) 1(^(2,0)^)1 = ^^. 

Proof. To construct if), we use the multiplicities in the spectrum: From 
the formulas ( |5.3|) for the eigenphases <f) V) i we see that <p v> i = 4> v ' t v mod 1 
if and only if (r/) 2 = rj 2 mod D, and in addition 

l' = l- V + b(r] - rj') mod M. 

In particular the multiplicity of <f> n j is exactly 

mod D' : (i]') 2 = rj 2, mod D}, 

which is independent of I. As a special case, we have 01,-6 = 0-i,& if 
D > 2. 
Now take 

(6.2) ^ = -^1,-6 + ^-1,6)- 
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Then since ipi-b, ip-i,b are orthonormal eigenfunctions with the same 
eigenphase, ip is a normalized eigenfunction. We compute 

(6.3) 

(T N {2,0)iP,iP) =^(7X2,0)^1,-6,^1,-6) + (^(2,0)^1,-6,^-1,6) 

+ (^(2,0)^-1,6,^1,-6) + ^(2,0)^-1,6,^-1,6) 



By Lemma |573| , we have (7^(2, 0)ip v ,i, Vv,'') = unless rf + rii = 
r) mod D. Thus in our case if D > 2 we see that all but the second 
summand in (16.31) are zero. As for that second summand, we see from 



Lemma |5.3j that in absolute value it equals 

1(^(2,0)^1,-6,-0-1,6)1 = 

where the Gauss sum G(—2b, —2b\ M) is given by 

-2b(x 2 + x) 



G(-2b, -2b; M) — e 



M 

xmodM 



In particular, if M is odd then by Lemma |4.5| its absolute value is v M. 
Thus we find that if M is odd then for our eigenfunction ip in ( |6.2| ) we 
have 

1(^(2,0)^)1 = ^, 
as required. □ 



6.2. Proof of Theorem We 

require the following construction: 



Lemma 6.4. Given an increasing positive function g(x), there is an 
irrational a, which has approximants \a — a/N\ < 1/N with arbitrarily 
large N so that M = N/ gcd(a, N) is odd, and satisfying 

M^>g{N) 2 . 



Once Lemma |67J] is proved, we will then take ip as in Proposition [673 



and then since M is odd, ( |6.1|) holds so that 

|(T^(2,0)V,V)| = 7T77T7?» ' 



2MV2 " g (N) ' 

which will conclude the proof of Theorem |6.2|. □ 
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6.3. Continued fractions. To prove Theorem |6.1| and Lemma f>.\ 
we first review some basic facts about continued fractions; see || for 
details. Give a sequence of integers ao, ai, 0,2, ■ ■ ■ with aj > 1 if i > 1, 
consider the (finite) continued fraction 

[ao; ai, 02, ... , a n ] = ao H ■ j — . 

CLi H 

The "partial convergents" p n , g n are defined through the recursion 
(6.4) p n+ i = a n+ ip n + p n -i, q n+ i = a n+1 q n + g n _ 1( n > 1 
with initial conditions po = a o, 9o — 1 5 Pi — a i a o + 1; Qi — a i- We have 

— = [ao; ai, 02, . . . , a n ]. 

The partial convergents satisfy the relation 

p n q n -i - Pn-iQn = (-l) n_1 , 

from which it follows that p n and q n are co-prime, and that q n -i and 
q n are co-prime. In particular at least one of q n -i, q n is odd. Another 
consequence is 

Pn Pn-l 



(6.5) 



Qn Qn-1 Qn-lQn 



We now construct the continued fraction a := [a ; ai, 02, . . . ] as the 
limit 

r 1 1 r Pn 

a := [ao; ai, 02, ... J = ao H 1 — = hm — 

ai H «^oo g n 

a2+ ' • 

(the limit exists by virtue of ( |6.5| )). It defines an irrational number. 

Conversely, for any irrational a, set ao = a, and for n > define 
integers a n and reals a n+ i > 1 by a n = a n + l/a n+1 . The integers a n 's 
are called the "partial quotients" of a, and are positive if n > 1. Then 

1 1 1 
a = [a ; ai, a 2 , . . . , a„_i H J = a H . 



ai + • + - — — T- 

a n-l + — 



We have 



and 



a = 

Oin+lQn + Qn-1 



Pn (-1)" 

a 



q n q n (a n+ iq n + q n - X ) 
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In particular, since a n+ i > a n+ i and q n > 1 one gets 
(6.6) |«_^|<^^. 



6.4. Proof of Theorem |6.1| . By Proposition |6.3|, given a it suffices 



to find arbitrarily large N and approximants \a — a/N\ < 1/N so that 
M = N/ gcd(a, N) is odd and satisfies N = M 2 . To do so, let ^ = J* 
be a partial convergent with M = q n odd. Since at least one of g n _i, 
q n is odd, there are infinitely many such M. Set 

N :=q 2 n = M 2 , a := p n q n , 

so that D := gcd(a, N) = q n . We have 



l a , , p n 1 1 



so all our requirements are satisfied. This proves Theorem IO. □ 



6.5. Proof of Lemma |6.4| . We begin with a construction of an irra- 
tional: 

Lemma 6.5. Given any positive increasing function F(x) there is an 
irrational a so that there are arbitrarily large q n and approximants 
Pn/q n so that F(q n ) < a n+x q 2 n and 

i ~Pn< 1 
\a < 



q n F{q n ) 
Moreover, we can require that q n are all odd. 

Proof. We define a through its continued fraction expansion, that is 
via the partial quotients a n . Set a = 0, and a± > 1 to be integer 
with ai > F(l). We define the partial quotients a; t inductively: Given 
do, di, . . . ,a n , we get the partial convergents p n , q n , and now choose 
a n+ i to be an integer so that a n+ i > F(q n )/q 2 . Set a = [a ;ai,a 2 , ■ ■ ■]■ 
Then from ( |6.6| ) 

\a < < 



q n a n+1 q 2 F(q n ) 

by our choice of a n+ \. Since at least one of every pair of consecutive g n 's 
is odd, we get infinitely many p n jq n satisfying our requirements. □ 



To conclude the proof of Lemma |6.4| , that is to find the required a, 
set G = g 2 , which is increasing. Then let F = G^ 1 be the inverse 
function to G which exists since G is increasing, and is positive. 
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Using Lemma |6.5| , we construct an irrational a whose partial con- 
vergent s satisfy 

i P» l / 1 ^ 1 

\a < - < — r. 

q n a n+1 qi F(q n ) 
Now take n so that q n is odd (there are infinitely many such n's) and 
set M := q n , b := p n (these are co-prime), and 

iV = a n+l ql, a = a n+ iq n p n 

so that D := gcd(a,N) = a n+1 g n , and \a — a/N\ < 1/N. Finally, 
M < G(N) = g{N) 2 because F(q n ) < a n+1 ql = N and since G is 
increasing, M = q n = G(F(q n )) < G(N) = g(N) 2 as required. □ 
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